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ABSTRACT 

The  equations  of  motion  for  steady  incompressible  flow  in  an  annulus 
with  fluid  removal  or  injection  at  the  porous  walls  have  been  reduced  to  a  pair 
of  ordinary  differentia*  equations  by  Terrill  (3) .  Solutions  for  the  case  in 
which  one  wall  rotates  are  not  yet  available;  however,  the  solution  of  the  rela¬ 
ted  problem  of  a  plane  channel  with  a  moving  wall  can  be  expressed  in  terms 
of  the  corresponding  results  (4,  5)  for  a  stationary  wall.  That  solution  is  ob¬ 
tained  in  the  present  work.  The  solution  is  used  to  determine  the  drag  on  the 
moving  wall.  Preliminary  experimental  results  with  an  apparatus  consisting 
of  a  porous  outer  cylinder  and  a  solid  rotating  shaft  sized  such  that  the  radius 
ratio  is  0  -926  indicate  satisfactory  agreement  for  the  case  of  fluid  injection; 
however,  there  is  poor  agreement  for  the  fluid  removal  case. 
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CHAPTER  I 


INTRODUCTION 


Phis  thesis  covers  an  analytical  study  of  the  drag  forces  acting  on  a 
shaft  rotating  concentrically  in  a  stationary  porous  cylinder,  through  which 


fluid  is  injected  or  withdrawn.  The  radius  ratio  of  the  annular  section  formed 
by  the  stationary  porous  cylinder  and  the  rota- 'tug  shaft  is  at  least  0.8  in  order 
that  the  annular  section  may  be  treated  as  a  plane  channel.  Design  of  an 
experimental  apparatus  for  testing  the  analytical  study  is  also  treated  and 
n  preliminary  experimental  results  are  presented. 

The  study  is  restricted  to  the  steady  laminar,  incompressible  flow  of 

: * 

a  Newtonian  fluid. 

Previous  work  will  be  cited  with  numbers  in  parentheses  which  refer 
to  the  List  of  References. 


A.  Previous  Work 


A  number  of  solutions  for  steady  laminar  flow  in  channels  with  porous 
walls  have  appeared  in  the  literature  over  the  past  few  years.  A.S.  Berman 
(i)  presented  the  first  solution  to  the  Navier-Stokcs  equations  for  steady 
laminar  flow  tn  a  uniformly  porous  channel,  Berman  (2)  also  gave  a  solution 
for  the  laminar  steady-ctate  Cow  of  an  incompressible  fluid  in  an  annulus 
formed  by  concentric  uniformly  porous  cylinders  for  the  special  case  where 
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the  rate  of  fluid  injection  through  ->ne  porous  cylinder  was  equal  to  the  rate  of 
fluid  withdrawal  at  the  other  porous  wall .  R.  M.  Terrill  (3)  reduced  the  pro¬ 
blem  of  laminar  flow  through  a  porous  annulus  with  constant  suction  velocity  at 
the  wall,  and  with  swirl  to  the  solution  of  four  nonlinear  differential  equations. 
Terrill  and  Shrestha  (4)  investigated  incompressible  laminar  flow  in  a  channel 
with  porous  walls  of  different  permeability  and  various  combinations  of  flow 
rates  through  the  wails.  Most  recently.  M.  S.  Tsai  (5)  obtained  solutions  to 
the  channel  problem  with  one  wall  porous  and  to  die  annulus  problem  with  one 
wail  porous.  Tsai  used  die  perturbation  method  (1J.  method  of  averages  (6j. 
and  numerical  techniques  to  obtain  solution  for  both  the  channel  and  annulus. 

He  compared  the  results  of  die  annulus  and  channel  solutions  and  used  the 
comparisons  to  determine  the  values  of  injection  and  suction  rates  and  the 
annulus  radius  ratio  for  which  the  annulus  may  be  treated  as  a  flat  channel. 

Tsai's  (a)  results  have  been  tested  experimentally  at  the  University  of 
Tennessee  by  Dr.  H.  Weissberg  and  Mr.  f.  Curlcc.  The  results  indicate  a 
good  correlation  between  Tsai’s  study  and  the  experimental  data. 

15.  Present  Work 

The  present  work  addresses  die  problem  of  extending  Tsai’s  to) 
investigation  to  the  case  of  a  rotating  shaft,  rather  than  a  stationary  one.  Use 
is  made  of  a  plane  channel  approximation  to  get  ..  novel  solution  of  the  Na*  icr- 
Stokes  equation  for  the  tangential  velocity  distribution  in  an  nnnutus  with  a 


CHAPTER  II 


EQUATIONS  OF  MOTION 

A.  Assumptions  and  Equations 

The  following  assumptions  are  made  in  treating  the  equations  of  motion 
for  a  shaft  rotating  within  a  porous  cylinder: 

1.  The  fluid  llow  is  laminar,  incompressible  and  fully  developed. 

3.  The  velocity  of  the  fluid  through  the  porous  wall  is  independent  of 
position. 

3.  The  flow  is  axisym metric;  the  shaft  and  porous  cylinder  are 
concentric. 

Tsai's  (5/  results  for  a  plane  channel  approximation  to  the  annulus 
problem  indicate  that  one  may  use  the  approximation,  with  less  than  5-percent 
error  in  the  calculated  pressure  drop,  provided  the  radius  ratio  of  the  annulus 
section  is  greater  than  0.9  and  the  wall  Reynolds  number  for  the  plane  channel 
is  less  than  2.0.  The  vvaii  Reynolds  number  is  defined  as  R  -  p  Yh/p ;  where  V 
is  the  radial  velocity  of  the  fluid  through  the  porous  wall  (positive  for  suction), 
h  is  the  distance  separating  the  plane  solid  wall  from  the  plane  porous  wall, 
u  is  the  viscosity  of  fluid,  and  p  is  the  mass  density  oi'  the  Ruin. 
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Tsai's  (5)  plane  channel  approximation  solution  for  the  annular  section 
was  with  a  stationary  solid  wall.  In  the  present  work  this  approximation  is 
generalized  to  treat  the  case  of  a  concentrically  rotating  solid  shaft.  The 
rotating  shaft,  as  shown  in  Figure  1,  corresponds  to  a  plane  solid  wall  moving 
in  a  direction  peroendicular  to  the  average  fluid  motion. 

A  set  of  Cartesian  coordinates  is  chosen  for  the  plane  channel  approxi¬ 
mation  such  that  the  x  coordinate  direction  corresponds  to  the  axial  direction 
of  the  annulus  section.  The  y  direction  corresponds  to  the  radial  direction  of 
th  ;  annulus.  The  moving  plane  solid  wall  is  physically  located  at  y  =  0  and 
corresponds  to  the  surface  of  the  shaft  in  the  annular  section.  The  z  direction 
completes  die  right-handed  orthogonal  coordinate  system  indicated  in  Figure 
1  and  corresponds  to  the  tangential  velocity  component  of  the  shaft  in  the 
annular  section. 

The  Xavier-Stokcs  equations  of  motion  for  the  steady,  laminar, 
incompressible  flow  in  the  plane  channel  arc 


/  3u  <*u  Su  \  3n  /  Ai  8ni  Oni  \ 


(  c3v  tV 

piur'-vr 
\  ox  &v 


ftv\  cV>  /  f rv  8V  fry  \ 

"te/*  *  +  |lU?  +  5**  a? ) 


V  9y  ?Z  /  Pz  *\ax2  dv2  dz2  ) 


where  u,  v,  and  w  are  the  fluid  velocity  components  in  the  x,  y,  and  z  coordi¬ 


nate  directions,  respectively,  and  p  is  the  fluid  pressure. 

The  equation  of  continuity  for  an  incompressible  fluid  is 


+  JK.  +  -xs  -  0 

ax  dv  fa 


The  axial  velocity  component  u  and  the  pressure  p  are  assumed  to  be 
functions  of  the  coordinates  x  and  y  only.  This  assumption  is  based  on  the 
shaft  rotating  concentrically  in  the  porous  cyi'nucr;  extending  this  to  the 
plane  channel  approximation  indicates  that  the  distance  h  separating  the 
stationary  porous  wall  and  the  solid  moving  wall  remains  constant  and  inde¬ 
pendent  of  coordinate  directions.  The  velocity  at  the  stationary  porous  wall 
is  assumed  to  l>e  perpendicular  to  the  wall  with  the  constant1  value  V  and  the 
velocity  component  v  is  assumed  to  he  a  function  of  the  y  coordinate  only. 

The  velocity  component  w  in  J:e  z  coordinate  direction  is  also  assumed 
to  be  a  function  of  the  coordinate  y  only.  The  assumptions  are  summarized  as 


follows: 


=  u  (x,  y). 


=  p  (x.  y) , 


=  v  (y). 
~  w  (y) . 


1.  In  practice  the  constant  velocity  through  the  wall  was  achieved  by 
insuring  that  the  pressure  drop  across  the  porous  wall  was  much  greater 
than  the  axial  pressure  diop  in  the  annulus. 
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Substituting  equation  (5)  into  equations  (1),  (2),  (3?,  and  (4)  leads, 
respectively,  to 


The  boundary'  conditions  are 


(X.  0)  -  0. 

V  (0)  =  ft,  w 

(01  =  W 

(x.  hj  =  o. 

v  (h)  ~  V,  and  \v 

(h)  =  ft  , 

(10) 

where  the  velocity  V  is  tire  injection  or  suction  velocity  at  the  porous  boundar 
and  W  is  the  velocity  of  the  moving  wall. 

Tsai  (5}  and  Terrill  (6j  reduced  equations  (6)  and  (7)  subject  to 
equation  to.)  to  the  dimensionless  ordinary  differential  equation 

F'"'  +  r(f’F!’  -  FF'”)  -•  0  ,  (il) 
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where  F  is  a  function  of  the  dimensionless  variable  X  =  y/h  and  R  is  the  wall 
Reynolds  number.  The  velocity  components  u  and  v  are  related  to  F  and  F'  by 

u  (x,  y)  =  (0)  -  y)  F’  (X) 

and 

v(y)  =  VF  (X)  ,  (12) 

where  U  (0)  is  the  average  velocity  in  the  x  direction  at  x  =  0.  The  boundary 
conditions  on  equation  (11)  are  found  by  using  equations  (10)  and  (12>. 


F(0)  = 

0. 

F (1 )  =  1 

F'(0)  = 

<>, 

F* (1)  =  0 

(13) 


Now,  in  addition  to  solutions  for  G  and  7.  the  solid  moving  wall 
requires  that  a  solution  be  found  for  equation  (8).  The  velocity  component  in 
the  z  coordinate  direction  is  assumed  to  be  of  the  form 


w  =  WG(X) 


(14) 


where  W  is  the  velocity  of  the  solid  wall  and  G(X)  is  a  dimensionless  function 
of  the  variable  X.  Substitution  of  equation  (14)  in  (Si  and  use  of  (12)  along 
with  Uie  definition  of  wall  Reynolds  number  R  =  yields 

p 


G"  (X)  =  R  F  (X)  G’  (X)  . 


(15) 
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The  boundary  conditions  on  equation  (15)  are  found  from  equation  (14) 
and  the  no  slip  condition  at  the  solid  and  porous  wails: 

G  (0)  -•  1,  G  (1)  =  0.  (16) 


The  solution  to  equation  (15)  subject  to  the  boundary  conditions  (16)  is 

X 


G(A)  --  1  - 


\  /  a  \ 

f  exp  I  f  RI*' (y)dy  Ida 
0  \  0  / 


1  /  0!  \ 

J  0KP  (  f  RF(y)dy  )da 
0  \  0  / 


(1?) 


It  is  interesting  that  die  reduction  of  equations  (6)  and  (7)  to  the 
dimensionless  ordhary  differential  equation  (11)  is  independent  of  the  motion 
of  the  wall  in  the  z  coordinate  direction;  hence,  the  resulting  velocity  compo¬ 
nents  in  the  x  ana  y  coordinate  dit\  tions  arc  valid,  regardless  of  the  motion 
of  the  'vail  in  the  z  coordinate  direction.  However,  equation  (S)  and  its  re¬ 
duction  to  equation  (15)  arc  dependent  on  the  solution  to  equation  (11).  It  has 
thus  been  shown  that  the  existing  solutions  of  equation  (11)  can  be  incorpora¬ 
ted  into  the  new  results,  equation  (17),  to  describe  the  case  of  shaft  rotation 
and,  indeed,  this  is  the  procedure  that  is  used  in  the  present  work. 


B.  Equation  for  Viscous  Drag 


The  shearing  stress,  i.  e-. ,  drag  force  in  the  z  direction  per  unit  area  of 


the  moving  wall,  is 


11 


Equations  (5)  and  (14)  are  substituted  into  equation  (18)  to  yield 


D  =  a  =  K, 


h 


KG’  (0) , 


(19) 


where  the  constant  K  is  defined  as 


W 


K  -  T  ' 


(20) 


The  expression  for  G’  (0)  is  found  from  equation  (17)  so  that  equation 
(19)  may  be  written  as 


D  = 


K 


J 


0 


f  <*P 


RF(y)  dy^da 


(21) 


It  is  interesting  to  take  a  closer  look  at  equation  (21) .  The  case  of  ao 
flow  through  the  porous  wall,  R  -  0,  corresponds  to  the  value  unity  for  the 
denominator;  hence,  the  ratio  of  the  drag  forces  without  flow  to  the  drag  forces 
with  flow  through  the  wall  is  just  the  value  of  the  denominator  of  equation  (21), 
that  is  l/cr  (0). 

The  viscous  drag  forces  acting  on  a  rotating  shaft  may  be  determined 
experimentally  by  measuring  the  time  rate  of  change  of  angular  velocity  pro¬ 
vided  that  all  external  forces  are  limited  to  those  associated  with  viscous 
drag.  An  expression  for  angular  momentum  decay  is 


L  -  -  J 


dt 


(22) 


where.  L  is  the  applied  torque  or.  in  this  case,  the  product  of  the  viscous 
force-  and  the  moment  arm,  J.  is  il*e  mass  moment  of  inertia  and  w  is  the 
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angular  velocity.  Equation  (19)  is  the  expression  for  the  viscous  drag  force 
per  unit  area  and  may  be  converted  to  an  expression  for  the  drag  torque  by  the 
appropriate  introduction  of  shaft  area  exposed  to  the  viscous  drag  and  the  mo¬ 
ment  arm;  hence, 

L  =  DAr  ,  (23) 

where  A  is  the  shaft  surface  area  and  r  is  the  shaft  radius.  Equation  (23) 
may  be  expanded  to 

L  =  (■'rH2"* 

However.  V*  is  the  velocity  if  a  point  on  the  shaft  and  is 

W  -  r.o  ,  (25) 

where  w  is  the  angular  rate  of  the  shaft.  Equatic  \s  (25)  and  (2A)  com¬ 
bined  become 


L  =  K{  Gf  (0)  uJ  .  (2G) 

where 


Since  equations  (22)  and  (26)  arc  both  expressions  for  the  shaft  torque  due  to 
viscous  forces,  they  may  be  combined  to  yield 


KjG'  (0)  -o 


(27) 


Equation  t27)  is  subject  to 


w  =  cu0  at  t  =  0  , 


where  w0  is  the  angular  rate  at  the  start  of  any  rundown. 

A  solution  to  equation  (27)  subject  to  the  boundary  condition  (28)  is 


*  u 
in  —  =  - 


In  terms  of  the  time  period  per  revolution,  r  =  2»/ u,  equation  (29)  may  be 
written  as 


Ki  C  (0) 


where  r0  is  the  initial  time  period  for  one  revolution  of  the  shaft  and  t  is  a 
continuous  variable  which  expresses  the  period  at  any  subsequent  time. 
Equation  (30)  provides  the  basis  for  the  experimental  work  to  be  discussed  in 
die  next  chapter. 


CHAPTER  III 


EXPERIMENTAL  WORK 


The  objective  of  the  experimental  part  of  the  work  is  to  make  a  preli¬ 
minary  investigation  into  the  validity  of  using  tne  previously  discussed  plane 
channel  approximation  in  solving  the  Navier-Stokcs  equations  for  a  solid  shaft 
rotating  in  a  porous  sleeve-  The  equipment  described  in  Tsai's  (5)  investiga¬ 
tion  was  available  for  this  work;  however,  modifications  were  necessary  to 
provide  for  rotation  of  the  shaft.  The  experimental  apparatus  was  modified 
accordingly  and  is  described  below. 

Our  experimental  test  of  the  present  thco»\  tical  work  consisted  of 
using  the  rundown  history  of  a  rotating  shaft  to  obtain  experimental  values  of 
G'  (0)  from  equation  (30) .  That  work  is  also  described  in  this  chapter. 


A.  Apparatus 


In  the  following  description  of  the  test  equipment,  reference  is  made  to 
the  diagram  of  the  apparatus  (Figure  2)  and  the  lettered  parts  designated  on 
the  diagram. 

The  ends  of  the  porous  sintered  titanium  cylinder  (E)  are  attached  to 
the  flanges  (F)  and  (G)  at  the  epoxy  joints  (J) .  These  flanges  are  attached  to 
flanges  (R)  and  (L),  thus  mounting  the  porous  cylinder  inside  the  jacket  HI], 


Figure  2 .  Diagram  of  apparatus 
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which  is  used  tc  supply  air  or  to  remove  it  from  the  outside  of  the  porous 
cylinder.  Three  rods  (i  ‘ ,  equally  spaced  around  the  flanges  are  used  to  pro¬ 
tect  the  thin  joints  (J)  from  rupture  in  case  of  exposure  to  excessive  air  pres¬ 
sure  in  the  jacket. 

The  rotating  shaft  (B)  is  supported  in  the  vertical  direction  by  the  gas 
thrust  bearing  (A).  Thrust  bearing  (A)  also  offers  some  lateral  support. 

The  uppermost  part  of  the  rotating  shaft  is  supported  radially  by  gas  bearing 
(C) .  Both  gas  bearings  (A)  and  (C .)  help  maintain  a  concentric  condition  of 
the  rotating  shaft  (B)  and  the  stationary  porous  cylinder  (E). 

The  annular  channel  consists  of  the  space  between  the  rotating  shaft 
(B)  and  the  porous  cylinder  (E) .  Calming  chambers  are  provided  by  the 
spacers  (K)  between  the  flanges  (M)  and  (R)  and  (M)  and  (L).  The  cham¬ 
bers  are  used  for  calming  the  air  entering  or  being  dischai'ged.  There  are 
thus  five  flow  paths  through  which  the  airstream  can  enter  or  leave  the  annular 
channel:  the  port  located  in  the  center  of  jacket  (fl).  the  two  ports  (not 
shown)  passing  through  the  spacer  rings  (K)  into  the  calming  chambers  at 
each  end.  and  the  gas  bearings  (A)  and  [C).  The  two  ports  labeled  (N)  are 
used  for  monitoring  the  pressure  in  the  calming  chambers. 

Rotating  shaft  (Bi  is  actuated  by  turbine  (D)  and  two  jets  not  shown. 
The  shaft  is  constructed  of  aluminum  and  is  machined  both  outside  diameter 
and  inside  diameter.  The  shaft  provides  a  radial  gap  of  o.OiiOS  in.  The 
radius  ratio  for  this  shaft  and  the  porous  cylinder  (E)  is  0.P2G3.  'Hie  shaft 
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(B)  wall  thickness  is  small  in  order  to  keep  tlie  inertia  to  a  small  value  and, 
hence,  makes  the  effects  of  the  viscous  drag  more  pronounced. 

The  bottom  gas  bearing  fed  by  (Q)  and  the  top  gas  bearing  fed  by  (P) 
were  designed  to  give  minimum  drag.  The  gas  bearing  drag  and  end  effects 
were  calculated  to  be  less  than  three  percent  of  the  maximum  drag  anticipated. 

Shaft  speed  is  determined  by  utilizing  a  dark  spot  cn  the  shaft,  a 
Mechanical  Technology,  Inc. ,  Fotonic  sensor  and  two  Hewlett-Packard  time 
interval  counters.  The  Fotonic  sensor  generates  one  pulse  per  revolution  of 
the  shaft.  The  two  Hewlett-Packard  counters  are  used  to  measure  the  time 
for  one  complete  revolution  of  the  shaft  and  the  time  between  shaft  period 
measurements.  The  shaft  rundown  history  is  obtained  by  recording  the  coun¬ 
ter  readings. 


B.  Description  of  Experiments 

Tests  were  run  with  fiitered  dry  air  as  the  test  fluid  for  wall  Reynolds 
numbers  ranging  from  +.'1.359  to  -5.241.  Mass  flow  rates  through  the  porous 
wall  were  adjusted  and  the  wall  Reynolds  numbers  were  calculated  from  these 
rates  and  the  laboratory  conditions.  An  example  of  these  calculations  is 
shown  in  the  Appendix.  The  calculated  wall  Reynolds  numbers  were  used  with 
Tsai’s  (5)  perturbation  method  to  provide  values  of  F  in  equations  <17) .  The 
results  were  used  to  determine  the  theoretical  value  of  the  denominator  of 
equation  (21)  and,  hence.  O’ (0)  for  the  flow  conditions  imposed  during  the 


experi  mcncs . 


'Hie  experiments  and  data  aequisition  were  designed  to  obtain  an  accu¬ 
rate  graph  of  equation  (30)  in  versus  time  for  wail  Reynolds  numbers  cor¬ 
responding  to  those  used  in  the  analytical  determination  of  G*  (0)  .2  Accurate 
measurements  of  shaft  rotational  period  and  elapsed  time  were  accomplished 
with  precise  electronic  counters  and  a  shaft  rotation  sensing  device,  as  des¬ 
cribed  previously.  The  mass  flow  rates  into  the  test  apparatus  u'ere  accounted 
for  with  calibrated  critical  orifices  and  electronic  mass  flow-  meters.  Fluid 
temperature,  room  temperature  and  atmospheric  pressure  were  monitored 
during  the  testing  and  were  used  in  the  calculations  for  mass  flow'  rates,  fluid 
viscosity  and  wall  Reynolds  numbers.  Tests  were  run  with  no  mass  flow' 
through  the  porous  wall,  mass  flow'  outward  and  mass  flow  inward.  These 
conditions  arc  described  on  the  following  pages. 

Positive  Reynolds  Numbers 

Positive  wall  Reynolds  numbers  correspond  to  the  case  w’hen  fluid  is 
allowed  to  flow  radially  outward  through  the  porous  cylinder  (E)  in  Figure  2 
(p.  15).  The  fluid  entered  the  system  through  bottom  bearing  and  was 
monitored  with  a  mass  flow  meter.  The  various  wail  Reynolds  numbers  were 
obtained  with  a  calibrated  critical  orifice  to  control  the  flow  into  the  bottom 
port  (N).  The  fluid  entered  the  annular  section  at  the  bottom  and  flowed 

The  equation  predicts  that  such  a  graph  should  be  linear  with  its 
slope  equal  to  -KjG’  (0)/-T, 
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axialiy  and  radially  outward  into  the  jacket  (H).  which  was  maintained  at 
atmospheric  pressure.  The  fluid  flowing  out  of  the  jacket  (H)  was  monitored 
with  a  mass  flow  meter  to  insure  that  only  the  fluid  flowing  into  the  annulus 
section  at  the  bottom  flowed  axially  and  radially  through  the  porous  cylinder 
■;  ).  Huid  flow  to  *  c  top  nearing  was  discharged  through  a  resistance  at  port 
(N)  and  was  adjusted  to  insure  that  no  fluid  from  the  top  bearing  entered  the 
annular  section  at  the  top. 


Negative  Reynolds  Numbers 

Negative  wall  Reynolds  numbers  correspond  to  fluid  flow  into  the 
jacket  tit.  ( Figure  2,  p.  15)  and  through  the  porous  cylinder  (E)  into  the 
annulus  formed  by  porous  cylinder  (E)  and  shaft  (B).  The  fluid  was  dis¬ 
charged  through  top  port  (N)  and  pons  in  the  tou  bearings  support.  Fluid 
for  the  bottom  bearing  was  withdrawn  from  the  bottom  calming  chamber  at 

bottom  port  (N)  and  was  monitored  with  a  mass  flow  meter  to  insure  that 
oniv  fluid  for  the  bottom  hearing  was  discharged  at  bottom  port  (N).  A  cali¬ 
brated  critical  orifice  controlled  the  flow  entering  the  jacket  (II). 


Zero  Reynolds  Number 

The  zero  wail  Reynolds  number  corresponds  to  zero  induced  fluid  flow 
in  the  annular  section  formed  by  shaft  (B)  and  porous  cylinder  (E)  in  Figure 
2.  Fluid  from  the  bottom  thrust  bearing  was  discharged  at  bottom  port  (N). 


Fluid  for  the  top  bearing  was  vented  through  ports  in  the  bearing  support 
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plate  and  also  through  top  part  (N) .  The  jacket  (H)  was  closed  to  insure  that 
no  fluid  flowed  through  the  porous  cylinder  (E) .  With  this  arrangement  there 
is  also  no  axial  flow.  The  pressure  inside  the  jacke'  (H)  and  calming  cham¬ 
bers  at  each  end  of  the  apparatus  was  monitored  to  insure  that  it  was  atmos¬ 
pheric  . 


The  test  procedure  for  ail  three  flow  conditions  was  as  follows: 

1 .  The  flow  conditions  were  adjusted  to  obtain  a  desired  w'all  Raynolds 
number.  Calibrated  critical  orifices  and  electronic  mass  flow  meters  along 
with  flow'  restrictors  were  utilized  to  accomplish  the  desired  flow'  condition 
for  each  experimental  run.  The  shaft  was  rotated  slowly  during  the  flow' 
adjustments  to  insure  a  concentric  relationship  between  the  shaft  (B;  and 
porous  cylinder  (E)  of  Figure  2  (p.  15).  The  fluid  temperature,  room 
temperature,  and  barometric  pressure  were  observed  and  recorded  for  each 
experimental  run. 

2.  The  shaft  was  rotated  and  allowed  to  stabilize.  All  controllable 
external  torques  were  removed  so  that  only  the  viscous  drag  imposed  by  die 
flow  conditions  remained.  The  shaft  speed  decay  measurements  were  initia¬ 
ted  using  one  Hewlett-Packard  counter  to  measure  the  time  period  of  one 
complete  revolution  of  the  shaft  and  another  to  measure  elapsed  time.  The 
successive  period  measurements  and  elapsed  time  were  made  using  the  gate 
closing  pulse  of  the  time  period  counter  to  actuate  die  second  counter  in  a 
time  interval  mode,  'flic  manually  opened  gate  of  the  time  period  counter 
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stopped  the  time  interval  counter  and  closing  of  the  gate  after  one  complete 
revolution  of  the  shaft  started  the  time  interval  counter  for  another  elapse 
time  measurement.  The  operation  of  the  counters  in  this  manner  allowed  the 
time  period  per  revolution  and  elapsed  time  between  period  measurement*  to 
be  determined  with  an  accuracy  of  100  user.  The  manual  selection  of  the  auto¬ 
matic  timing  sequence  allowed  the  operator  ample  time  to  record  the  shaft 
period  and  elapsed  time  before  starting  the  next  sequence.  This  procedure 
resulted  in  10  to  20  sec  between  shaft  period  measurements.  The  total  time 
elapsed  for  any  experimental  run  was  from  200  to  300  sec.  The  shaft  speed  at 
the  beginning  of  the  test-  runs  ranged  from  about  200  to  300  rpm. 

3 .  Each  flow  condition  was  tested  with  the  shaft  rotating  first  in  the 
clockwise  direction  and  then  in  the  counterclockwise  direction.  This  insured 
that  any  asymmetry  associated  with  rotational  direction  would  be  duly  noted. 

C.  Experimental.  Data  Treatment 

The  objective  in  the  treatment  of  the  raw  data  is  to  generate  an  accu¬ 
rate  shaft  speed  variation  with  time  such  that  an  experimental  graph  of  equa¬ 
tion  (30)  can  be  obtained  for  each  flow’  condition. 

The  output  of  the  shaft  rotation  f  easing  device  is  typically  a  square 
wave  electrical  signal  as  indicated  in  Figure  3.  The  raw  data  from  the 
experimental  runs  are  the  r's  and  T's  indicated  in  the  figure.  The  determi¬ 
nation  of  ti  e  t's  also  indicated  in  the  figure,  is  necessary  in  order  to  get  an 
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Figure  3.  Fotonic  sensor  waveform  and  timing  scheme. 


accurate  graph  of  the  shaft  speed  variation  with  time.  Figure  3  shows  that 


t 

n 


(31) 


A  small  computer  was  used  to  obtain  values  of  t  from  equation  (31) 
and  to  compute  die  corresponding  left  side  of  equation  (30),  The  determina¬ 
tion  of  die  mass  moment  of  inertia  J.  which  appears  in  equation  (30).  is  dis¬ 
cus  sod  in  the  Appendix. 

A  regression  analysis  was  used  to  fit  a  straight  line  to  the  datafn  t0/t 
versus  time  for  each  experimental  run.  This  procedure  gives  the  slope  and 
intercept  of  the  fitted  line  and  also  the  coefficient  of  a  0:1 -percent  confidence 
interval  for  these  two  statistical  parameters. 
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The  t-distribution  was  used  for  the  95 -percent  confidence  interval 
determination.  This  was  based  on  the  generally  accepted  value  of  30  data 
points  above  which  the  normal  distribution  is  used  and  below  which  the 
t -distribution  is  used.  Most  of  the  experimental  runs  in  the  present  work  con¬ 
tained  less  than  30  data  points . 

D.  Results 

The  experimental  values  for  G’  (0)  were  obtained  from  equation  (30) 
using  the  experimental  slopes  to  obtain  values  of  -KjG'  (0)/J.  ITte  slopes 
were  deter;  .ined  from  a  least  squares  straight  line  fit  of  the  experimental  data 
as  previously  explained.  The  right  side  f  equation  (30)  represents  the  slope 
and  the  time  variable  t;  hence,  the  coefficient  of  t  was  set  equal  to  the  slope 
and  the  drag  ratio,  G’(0),  was  calculated  for  each  experimental  run.  A  typi¬ 
cal  example  of  this  calculation  is  shown  in  the  Appendix. 

A  total  of  30  experimental  runs  were  made  at  nine  different  wall 
Reynolds  numbers.  Table  I  is  a  comparison  of  the  theoretical  and  experi¬ 
mental  results  for  the  30  runs.  Repeatability  of  the  data  was  good  within  5 
percent  for  more  titan  one  run  at  the  same  wall  Reynolds  number.  Several 
different  tans  were  made  on  different  days  at  wall  Reynolds  numbers  equal  to 
zero  and  Table  I  indicates  that  the  repeatability  is  excellent. 

Figures  4  through  12  illustrate  typical  runs  for  various  Reynolds 
numbers.  Each  figure  shows  the  theoretical  line  calculated  from  equation 
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Figure  4.  Shaft  speed  decay,  R  =  0. 
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Figure  5.  Shaft  speed  decay.  R  =  -0.9248. 
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Figure  6.  Shaft  speed  decay.  R  =  -1.8496. 
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Figure  7.  Shaft  speed  decay.  ft=  -3.0S27. 


RUN  NO.  ! 

REYNOLDS  NO.! 
SLOPE 
INTERCEPT 
MAX.  rpm  ! 


CLOCKWISE 

53-1 _ ‘ 

-5.2407 


COUNTER- 

CLOCKWISE 

53-3 

-5.2407 


ggitiMTl 


-0.0220+0.0059 


222 


CLOCKWISE 

COUNTERCLOCKWISE 

THEORETICAL 


40  €0  80  100  120  140  160  ISO  200  220  240  260  280  300  320  340  360 

TIME  (sec) 


Figure  9.  Shaft  speed  decay,  K  =  -5. 2407 
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figure  11.  Shaft  speed  decay,  R  -  +2. 1255. 
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(30)  and  the  least  squares  line  as  fitted  through  the  plotted  experimental 
points.  There  are  two  such  experimental  lines  because  data  were  taken  for 
two  directions  of  rotation  as  explained  in  a  previous  section.  Description  of 
Experiments.  A  table  inserted  in  each  figure  also  provides  the  following  in¬ 
formation: 

1.  The  wall  Reynolds  number  for  the  run. 

2.  The  least  squares  intercept  and  its  deviation  (95-  percent  confi¬ 
dence  inte  rval)  for  each  direction  of  shaft  rotation. 

3.  The  least  squares  slope  and  its  deviation  (95-percent  confidence 
interval)  ior  each  direction  of  shaft  rotation. 

4.  The  shaft  speed  at  the  beginning  of  data  acquisition. 

The  results  of  the  experimental  runs  are  conveniently  summarized  by 
plotting  the  theoretical  (see  section  on  Viscous  Drag)  and  experimental  values 
(Table  I,  p.  24)  of  the  drag  ratio  GY  (o),  versus  wall  Reynolds  numbers;  such 
a  plot  is  shown  in  Figure  13.  The  computed  values  of  G’  (0)  are  given  for 
both  clockwise  and  counterclockwise  shaft  rotation  directions.  Reasonably 
good  agreement'1  between  the  theoretical  and  experimental  values  of  G’  (0)  is 
indicated  for  clockwise  shaft  rotation  with  wall  Reynolds  numbers  from  -1  to 
-5;  the  range  of  reasonable  agreement  for  counterclockwise  shaft  rotation 
extends  from  -1  to  -3.  Examination  of  Figures  5,  6,  and  T  (pp.  25-26)  indi¬ 
cates  that  the  slope  for  both  rotating  directions  in  the  wall  Reynolds  number 

3.  A  statistical  test  of  the  experimental  data  is  discussed  in  the 
appendix. 
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range  -1  to  -3  are  within  10  percent  oi  each  other.  Positive  wall  Reynolds 
numbers  resulted  in  no  apparent  agreement  between  the  theoretical  and 
experimental  results.  This  lack  of  agreement  >s  quite  evident  in  Figure  13. 
A  discussion  of  this  is  given  below. 


E .  Discussion 

The  lack  of  agreement  between  the  theoretical  predictions  and  the 
experimental  results  noted  aboye  suggest  that  an  instability  may  exist  for 
positive  wall  Rejaioids  numbers.  During  the  experimental  runs  with  positive 
wall  Reynolds  numbers,  it  was  observed  that  a  slight  induced  vibration  of  the 

apparatus  would  cause  the  shaft  to  liecome  unstable*  and  to  stop  (because  of 
contact  between  the  shaft  and  eidicr  of  die  bearing  stops)  in  a  few  revolutions. 
This  is  thought  to  be  associated  with  the  stiffness  of  the  shaft  air  bearings.  In 
order  to  obtain  any  d  un  over  a  reasonable  shaft  rundown  time,  it  was  neces¬ 
sary  to  carefully  accelerate  die  shaft  to  die  starting  speed  of  from  20ft  to  300 
m  and  see  if  die  rotation  was  stable  before  data  were  taken.  This  proce¬ 
dure  was  successful  as  the  slope  confidence  interval  shows  (Table  I,  p.  24). 
Although  the  experimental  work  was  preliminary  in  nature,  it  appears  unlikely 
that  the  large  uiscrcpancics  indicated  in  Figure  13  for  positive  wall  Reynolds 
numbers  were  introduced  by  die  experimental  apparatus.  It  is  speculated  that 

4.  Similar  instability  was  also  noted  for  negative  Reynolds  numbers  at 
starting  shaft  speeds  above  about  300  rpm. 


33 


an  instability  phenomenon  similar  to  Taylor  vortices  exists.  A  qualitative 
justification  for  this  speculation  follows.  The  critical  Taylor  number  (6)  is  an 


indication  of  vortex  formation: 


Wh  JT 


wh  re  W  is  the  peripheral  shaft  speed,  h  is  the  annulus  gap,  r  is  the  rotating 
shaft  radius,  and  v  is  the  kinematic  viscosity  of  the  fluid. 

In  the  present  work  the  maximum  shaft  peripheral  speed  was  19.1  in. 
per  second.  This  shaft  speed  resulted  in  a  critical  Taylor  number  of  8. 38  and 
would  indicate  a  stable  condition  at  least  for  R  =  0.  The  critical  Taylor 
number  could  not  be  calculated  for  positive  or  negative  wall  Reynolds  number 
since  Taylor  assumed  the  radial  velocity  component  to  be  zero;  however,  the 
circumferential  velocity  distribution  in  the  radial  direction  indicates  that  the 
formation  of  Taylor  vortices  ;s  plausible.  Figure  14  is  a  graph  of  the  cir¬ 
cumferential  velocity  distribution  m  the  radial  direction  for  both  positive  and 
negative  wall  Reynolds  numbers.  The  velocity  distribution  for  positive  wall 
Reynolds  numbers  appears  similar  (for  over  50  percent  of  die  gap)  to  what 
could  be  expected  for  concentric  rotating  cylinders  with  a  larger  radial 
clearance  that  was  used.  The  critical  Taylor  number  varies  as  the  gap  to  the 
3/2-  power  of  the  radial  clearance  according  to  equation  (32).  These  consid¬ 
erations  indicate  that  the  critical  Taylor  number  quite  possibly  was  reached 


for  positive  wall  Reynolds  numbers.  On  the  other  hand,  the  same  kind  of 
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Figure  1*J.  Velocity  distribution  in  annulus 
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speculation  would  predict  stability  for  negative  wall  Reynolds  numbers  since 
the  velocity  distribution  appears  similar  to  what  would  be  expected  with  a 


much  smaller  gap. 


CHAPTER  IV 


GG 

SUMMARY  AND  RECOMMENDATIONS 

A.  Summary 

The  present  work  addressed  the  problem  of  using  a  plane  channel 
approximation  to  obtain  theoretical  evpiessions  for  the  drag  forces  associated 
with  a  shaft  rotating  in  a  porous  cylinder  with  fluid  injection,  negative 
Reynolds  numbers,  cr  fluid  withdrawal,  positive  Reynolds  numbers.  The 
solution  to  the  Navier -Stokes  equations  for  fluid  Row  in  the  axial  and  radial 
direction  was  shown  to  be  independent  of  the  rotational  motion  of  the  shaft 
when  the  plane  channel  approximation  is  used:  it  predicted  the  linear  rela~ 
tionship  of  f  n  with  time  in  ’’coastdowr  .  xperiments. 

T 

A  statistical  analysis  of  the  preliminary  test  results  confirmed  the 
above  linearity.  The  experimental  data  agreed  reasonably  well  with  the 
theoretical  predictions  for  negative  wall  Reynold's  numbers:  however,  the 
experimental  data  for  positive  wall  Reynolds  numbers  did  not  agree  with  the 
theoretical  predictions.  It  was  speculated  that  the  lack  of  agreement  between 
theoretical  and  experimental  results  for  positive  wall  Reynolds  numbers  was 
due  to  an  instability  similar  to  the  formation  of  Taylor  vortices.  However, 
the  stiffness  properties  of  the  shaft  support  air  bearing  (and  associated  shaft 
instability)  seemed  to  be  such  that  any  disturbance  (c.  g.,  hydrodynamic 
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instability)  external  to  the  shaft  would  preclude  stable  rotation  of  the  shaft. 
This  possibility  is  considered  in  the  following  section  on  recommendations  for 
future  work. 


B.  Recommendations 

The  experience  gained  during  this  work  and  the  experimental  results 
suggest  a  continuation  of  this  work  in  the  positive  wall  Reynolds  number 
region.  A  lack  of  agreement  between  theory  and  experiment  was  shown  for 
positive  Reynolds  numbers;  the  shaft  support  design  was  not  suitable  for 
further  investigation  of  possible  causes,  e.  g..  the  previously  mentioned 
speculation  on  Taylor  vortex  formation.  The  use  of  precision  instrument 
bearings  cxibising  low  constant  drag  would  result  in  a  stiff  shaft  suspension 
system  that  would  allow  operating  speed  in  excess  of  1000  rpm  (much  higher 
than  our  present  300  rpm  maximum) .  Such  higher  speeds  would  allow  an 
invcstigatii  of  the  drag  in  the  vicinity  of  the  critical  Taylor  number.  Addi¬ 
tionally.  precision  instrument  bearings  would  insure  that  an  r.xisymmefcric 
condition  exists  between  the  shaft  and  the  porous  cylinder.  A  small  gap  h 
would  delay  the  onset  of  Taylor  vortices;  however,  a  porous  cylinder  with  a 
more  uniform  diameter  than  tee  one  used  in  this  work  (0.003  in.  out  of 
round  with  a  0. 050-in.  annular  gap)  would  be  necessary. 

The  plane  channel  approximation  applied  to  the  case  of  a  shaft  rotating 
concentrically  in  a  porous  cylinder  predicted  a  lessening  of  the  drag  on  the 


shaft  when  fluid  is  withdrawn  through  the  porous  cylinder.  It  would  be 
interesting  to  investigate  further  whether  there  are.  indeed,  any  conditions 
for  which  this  theoretical  prediction  could  be  confirmed  by  experiment  to  the 
same  extent  as  the  present  work  has  confirmed  the  predictions  for  fluid 


injection. 
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APPENDIX 


A.  Shaft  Moment  of  Inertia 


Equation  (30)  indicates  that  the  moment  of  inertia  should  be  known 
precisely  in  order  that  G'  (0)  be  determined  accurately  from  the  experimental 
value.  The  shaft  inertia  was  determined  with  the  aid  of  two  methods. 

1 .  The  shaft  dimensions  were  held  rigidly  during  machining  and  the 
final  dimensions  were  known  precisely.  This  information,  along  with  a  weight 
measurement  of  the  shaft,  was  used  to  calculate  the  mass  moment  of  inertia 
of  the  shaft. 

2.  The  shaft  moment  of  inertia  was  determined  experimentally  with 

the  aid  of  a  method  by  H.  M.  Angel  (7) .  The  shaft  was  suspended  by  three 

strings  and  caused  to  oscillate  about  its  longitudinal  axis.  The  weight  of  the 

shaft,  the  suspension  string  length,  the  shaft  radius  and  the  period  of  ocilla- 

tion  were  used  to  calculate  the  moment  of  inertia  of  the  shaft.  The  equation 
WR^  P“ 

used  is  -J  -  .  v; —  ,  where  J  is  the  mass  moment  of  inertia  of  the  shaft  in 
4tT*  L* 

in.  -lbj.-sec2.  W  is  the  shaft  weight  in  pounds,  P  is  the  oscillation  period  ir. 
seconds,  and  It  and  L  are  the  suspension  radius  and  length,  respectively.  A 
typical  trial  resulted  in  these  data: 


=  0.5105  Ibj. 

=  0.7320  in 
=  0.5226  sec 
=  4.5  in. 
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and  a  J  of  4.209  x  io“4  in.-lb^-sec2.  A  small  ring  used  to  support  the  shaft 
had  a  calculated  value  of  0. 13ft  x  10"*1  in. -lb^-s^c2.  The  resulting  shaft  mass 
moment  of  inertia  was  4  .074  x  10~4  in.  -Ib^-sce2.  The  calculated  value  as 
determined  by  item  1.  above  was  4. 130  x  io“4  in.  -Ib^-sec2  or  a  difference  of 
0.056  x  10  4.  This  represents  a  possible  error  in  measurement  of  less  than 
2  percent. 


B.  Sample  Calculations 


i.  The  wail  Reynolds  number  can  be  expressed  by  the  following: 


R 


A  typical  calculation  using  the  condition  imposed  on  run  60-1  is 

h  =  0.050S  in. 

m  -  0.0036317  lb  /sec 
m 

u  -  1.0133  n  io-G  lb  in. -sec 
m 

A  -  54.20S  in.2 
R  =  3.35S7  . 

2.  The  values  for  O’  (0)  were  calculated  from  the  experimental  data. 
A  typical  calculation  using  the  data  from  run  60-1  is 

(6)  ~  (experimental  slope)  .  1.5656 

Experimental  slope  -  -0.0039/sec. 


Kt  =  .1.01475  x  lo“c  in-ibf-sec 
J  =  4.074  x  10“4  in-lb^-sec2  . 


3.  The  critical  Taylor  number  may  be  expressed  by  the  following: 


_  wh  nr  , 

T  =  -  -  =  41.3  . 

a  v  r 


Run  55. 1  data  may  be  used  for  this  calculation  and  are 

W  *  19.1  in. /sec 


r  =  0.C386  in. 


h  =  0.0508 


v  =  0.03267  in. 2/scc  . 


The  resulting  Taylor  number  with  these  data  is  8.376. 


C.  Statistical  "Luck-of- Fit-Test*' 


An  F-test  was  conducted  on  the  data  to  test  the  hypothesis  that  the 
experimental  results  were  in  good  agreement  with  the  theoretical  predictions 
for  G*  (0).  This  test  used  the  ratio  of  the  sample  variance  of  the  theoretical 
and  experimental  values  of  G’  {o)  aud  the  properly  weighed  sample  variance 
of  the  snaft  rundown  data.  The  results  of  the  test  indicated  that  the  differ¬ 
ences  between  the  theoretical  and  experimental  values  of  G1  (o)  (Figure  13, 
p.  31)  for  our  ranges  of  ”r  sonable  agreement”  are  not  due  to  random 
experimental  errors  in  the  rundown  data.  Sources  of  experimental  bias  in 


the  apparatus  (which  was  not  originally  designed  for  shaft  rotation)  were, 
course,  expected.  The  magnitude  of  such  bias  is  perhaps  indicated  by  the 
difference  between  the  results  for  clockwise  and  counterclockwise  shaft 
rotation. 


